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SUMMARY 

The heats of fusion of several fluoride salt mixtures with melting points 
greater than 973 K have been estimated from a coupled analysis of the availa- 
ble thermodynamic data and phase diagrams. Simple binary eutectic systems 
with and without terminal solid solutions, binary eutectics with congruent 
melting Intermediate phases, and ternary eutectic systems were considered. 
Several combinations of salts have been Identified, most notably the eutectics 
L1F-22CaF2 and NaF-60MgF2 which melt at 1039 and 1273 K respectively which 
possess relatively high heats of fuslon/gm (>0.7 kJ/g). Such systems would 
seemingly be Ideal candidates for the light weight, high energy storage media 
required by the thermal energy storage unit In advanced solar dynamic power 
systems envisioned for the future space missions. 


INTRODUCTION 

Thermal energy storage (TES) subsystems will be required for advanced 
solar dynamic power systems to provide electric power for future space mis- 
sions during periods of low or nonexistent solar Insolation. While batteries 
or fuel cells can supply the necessary energy, they Impose a large weight 
penalty. Less heavy alternatives have been suggested, and these generally 
Involve the use of a material to store heat for an engine turning a genera- 
tor. For example, solar radiation can be utilized to melt a solid while the 
solar dynamic power system Is In sun light, and then this energy can be 
extracted from the solidifying melt when the power system Is In eclipse 
(refs. 1 and 2). 

One of the prime criteria for selection of a solid to liquid phase change 
TES material Is a high heat of fusion. Fluoride salts of Na, LI, Mg, K, and 
Ca have the highest enthalpies of fusion among all the molten salts; however 
they can only be used If the operating temperature of the heat engine Is com- 
patible with their melting points. This limitation can be removed If eutectic 
mixtures of two or more fluoride salts are considered. Unfortunately, data 
are lacking on the heats of fusion for most eutectics; therefore, as a first 
point for the selection of potential salt systems, the heats of fusion must be 
estimated . 

This report describes the estimation of latent heat of fusion of several 
multicomponent fluoride eutectics as determined by a coupled analysis of the 
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available thermodynamic data and the phase diagrams. All systems chosen have 
melting points above 973 K to be consistent with the potential operating tem- 
peratures of Brayton and Sterling cycle engines proposed for the advanced 
solar dynamic power systems. 
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SYMBOL LIST 

constants In the expression (eqs. (4), (5), and (7)) 
activity of compound 1 In the solution 

heat capacity of a solid compound 1, JK -1 mol -1 

heat capacity of a compound 1 In the liquid state, JK -1 
mol -1 

eutectic 

excess partial molar Gibb's free energy for compound 1 In 
the solution, kJ/mol 

charge of the cation 

excess entropy of mixing In the solution, JK _1 mol^ 1 
solid 

dry temperature, K 
melting point, K 

mole fraction of a compound 1 In a solution 

equivalent Ionic fraction of a compound 1 In a solution 

mole fraction of compound 1 In the a solid solution that 
Is In equilibrium with the eutectic 

mole fraction of compound 1 In the 0 solid solution that 
Is In equilibrium with the eutectic 

activity coefficient of compound 1 In the solution 

activity coefficient of compound 1 ' In the eutectic melt 

activity coefficient of compound 1 In the a and 0 solid 
solution, respectively 

Gibb's free energy of fusion of a compound 1 at any 
temperature other than the melting point, kJ/mol 
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heat of fusion of the eutectic, kJ/mol 

heat of mixing In the melt at the eutectic, kO/mol 

heat of mixing In the binary A-B and ternary A-B-C melt 

at the eutectic composition, kJ/ mol 

heat of fusion of a compound 1 at the melting point, kJ/mol 

partial molar enthalpy of mixing for the compound 1 In the 
solution, kJ/mol 

heat of fusion of a compound 1 at any temperature other than 
the melting point, kJ/mol 

heat of mixing In a mole of solution, kJ/mol 

heat of mixing per equivalent of solution, kJ/equIvalent. 

heat of mixing In the o solid solution, kO/mol 

heat of mixing In the B solid solution, kJ/mol 

heat of formation of the Intermediate phase e In the binary 
system A-B from pure compounds A and B, kJ/mol 

entropy of fusion of the eutectic, JK - ! mol -1 

Ideal entropy of mixing for the eutectic liquid, JK - ! 
mol -1 


entropy of fusion of compound 1, JK -1 mol -1 
total entropy of mixing, JK - ! mol - ! 

Ideal entropy of mixing In a solution, JK -1 mol - ! 

Ideal entropy of mixing for a solid solution that Is In 
equilibrium with the eutectic, JK -1 mol - ! 

Ideal entropy of mixing for B solid solution that Is In 
equilibrium with the eutectic, JK -1 mol - ! 

Intermediate compound 
liquid 
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THERMODYNAMIC CONSIDERATIONS 


General Approach 


The heat of fusion of any eutectic, AH e f, can be directly calculated via 
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or If the entropy of fusion, AS e f, can be estimated for 

the eutectic, 

where 
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4H ef * ss ef 

xT e 


(2a) 


where T e Is the eutectic temperature. Unfortunately, most of the terms on 
the right hand side of these equations have not been experimentally measured; 
hence they must be estimated through standard thermodynamic models and expres- 
sions. In this study AH e f Is computed by both equations (1) and (2a), and 
a comparison made; close agreement of the two values would Indicate a reliable 
estimate. 


Calculation of Terms for Equations (1) and (2) 

Heats of fusion of pure materials are only known at their melting points, 
but they can also be calculated for other temperatures T from 


<1 ’ ♦/ ( C M - c p.l) 


dT 


(3) 


T M 

where A H f , and AH, , are the heats of fusion of component 1 at tempera- 

' » 1 T * ' o S 

ture T and the melting point T m respectively; and Cp ^ and Cp ^ are 

0 

the heat capacities of the liquid and solid phases. While C , generally 

S P » * 

varies little with temperature, C p ^ is a strong function of temperature 
(ref. 3) 
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where A', B 1 , and C' are constants available from standard tables (ref. 4). 

The heat of mixing AH^ of many binary fluoride melts have been deter- 
mined by Kleppa and co-workers (refs. 5 to 8) and fitted to a polynomial 
expression of the type 

iH H “ ¥b ( A ' - B ' X B * C ' X B * D ' X b) < 5) 

where X A and Xg are the mole fractions of substances A and B, respec- 
tively In the mixture and A' f B' f C 1 , and D 1 are the constants. 

The entropy of mixing AS M for a binary melt Is the sum of the Ideal 

i j n r 

entropy of mixing AS^ U plus the excess entropy of mixing sjjj where 

4s’ d . -R[ X A l" X A ♦ X B m X B ] (6) 

and as suggested In reference 9 

S M ' X A X 8 (*' f B ' X B * CX B * °' X b) (7) 

with R being the universal gas constant. In general the excess entropy of 
mixing In the melt will only be Important If the enthalpy of mixing Is large; 

otherwise sjjj can be assumed to be zero. 

Equations (5) to (7) can also be used to describe the thermodynamic prop- 
erties of the solid solutions In which all the cations have the same valence. 
However, In charge asymetrlc Ionic solid solutions, the equivalent ionic frac- 
tion Instead of the simple mole fraction is used to calculate the enthalpy and 
entropy of mixing (ref. 10). The equivalent Ionic fraction X^ of a cation 
1 Is 
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where q^ Is the absolute charge of the cation 1 and the summation Is con- 
ducted over all the cations In the solution. Substitution of the equivalent 
Ionic fractions for A and B In equations (5) to (7) gives the enthalpy and 
entropy of mixing per equivalent of the binary solution. 

When thermodynamic data for the solutions are not available, frequently 
estimates based on the phase diagrams (refs. 9 and 11) are made by first cal- 
culating the Gibbs free energy of fusion of a component at a given temperature 
along the llquldus and then relating this to the excess Gibbs energy of the 
component In solution via 
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(9) 


RT Sin X* + G^> - RT in - G^ (S) = -AG^J* 

where G^ 8 ^ Is the excess partial molar Gibb's free energy for 1, AG°^ Is 

the free energy of fusion 1 at the temperature corresponding to the llquldus 
point. AG° can be determined from 


AG° ( P = aH T 
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f,1 
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The partial molar excess Gibb's free energy Is given by 


G^ = AH 1 - TS^ 


( 10 ) 


(ID 


where AH^ and are the partial molar enthalpy of mixing and excess par- 

tial molar entropy of mixing, respectively, for component 1. 

Each term In equation (11) can be related to simpler thermodynamic con- 
cepts via 

a 1 = X 1 • Yi (12) 

and 

G^ = RT an Yjj (12a) 


where a^ and Y ^ are the activity and activity coefficient of 1, respec- 
tively; and for binaries the application of Glbbs-Duhem equation gives 


dAH M 

AH 1 = AH H + (1 - X,) 


(13) 


S H * 
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For ternary systems the heat of mixing can be estimated from the equation 
derived by Olson and Toop (ref. 12) which relates the ternary heat of mixing 

A-B-C 

for a liquid solution consisting of A, B, and C, AH H , as a function of 

the binary values of AH U along composition paths with constant X./X D , X D /X_ 

rl A D D L 

and X./X_ ratios: 

A C 
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where the meanings of the terms shown In equation (15) are graphically shown 
In figure 1 . 


APPLICATION OF THERMODYNAMIC PRINCIPLES TO THE ESTIMATION OF HEAT OF FUSION 

In the following determinations the entropy of fusion of a component and 
the constants In the equations for enthalpies and entropies of mixing will be 
assumed to be temperature Independent. For simplicity four categories of 
eutectics are considered: binary systems (1) with and (2) without terminal 

solid solutions or (3) with congruently melting Intermediate compounds, and 
(4) ternary systems having no solid solutions. 


Binary Eutectic Systems Without Terminal Solid Solutions 


Heat of mixing known for the liquid . - The heat of fusion for the eutectic 

AF„ -BF„ (AH^~ B ) can be directly calculated from the enthalpies of fusion for 
q A q B V ef / / T T \ 

each salt at the eutectic temperature (fiH^ e A , AH^gland the heat of mixing of 

the liquid at the eutectic composition (aH^^ via 
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(16) 


Heat of fusion from entropy considerations . - The change In entropy during 
the melting of a eutectic can be expressed by 


AS*: 6 = X® AS* . + X® AS* „ + AS qU 
ef A f,A B f,B eM 


1d(S.) + $ E(l) 


eM 


(17) 


Available heat of mixing data for the fluoride melt systems that do not exhibit 
an Intermediate compound. Indicate that only small deviations from Ideality 
occur; hence the excess entropy of mixing will be assumed to be negligible. 

With Ideal entropy of mixing calculated from equation (6), the heat of fusion 
follows directly from equation (2a). 

Heat of mixing from phase diagram analysis . - Assuming excess entropy of 
mixing to be zero, the enthalpy of mixing of the salts at the eutectic can be 
calculated from the following equations: 
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where AG, A and AG, D are the Gibb's free energy of fusion of AF 
r,A r,b q. 

6 G 

and BFq , respectively at the eutectic temperature; y a and y b are the 
activity coefficient of AF and BF_ , respectively In the liquid at the 


<B 


,AB 


eutectic composition. Substitution of AH gM In equation (16) gives the heat 
of fusion of the eutectic. 

Example calculations . - Kleppa and co-workers have measured the heats of 
mixing for three systems (refs. 5 and 6) that exhibit a simple eutectic: 
KF-NaF, L1F-CaF2 and NaF-CaF 2 where 


AHj; F - NaF . -334 X RF • X Nap J/mol 


LIF-CaF, 


AH 


AH 


M 

NaF-CaF, 

< 

M 


= X L1F * X CaF 2 (- 5271 - 376 X CaF 2 


' X NaF * X CaF 2 (- 19413 " 376 X CaF 2 


J/mol 


J/mol 


( 21 ) 

( 22 ) 

(23) 


Estimated heats of fusion utilizing each of the above strategies for the sim- 
ple binary eutectic systems are given In table II; the close agreement among 
the values Is taken to signify a valid result. 


Binary Eutectic Systems With Terminal Solid Solutions 

Theory . - Blrchenall and Rlechman (ref. 13) have described a technique to 
estimate the enthalpy of fusion of the eutectic In binary systems with termi- 
nal solid solutions such as that shown In figure 2. At the eutectic tempera- 
ture, where equilibrium Is maintained among liquid of composition X e and 
solid solutions of composition X a and Xg, the heat of fusion per mole of 
eutectic Is 



(1 


X ) AH, e + X AH f e + AH^ U) 
e f,A e f,B eM 




(24) 


AB 

where AH^ Is the heat of mixing of the a solid solution at composition X 
AB 

and aH bh Is the heat of mixing of the B solid solution at composition X^ . 
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Similarly, assuming the excess entropy of mixing to be negligible, the entropy 
of fusion per mole of the eutectic can be written as 


AS 


AB 

ef 
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V AS f,A + X e AS f,B + AS eM 
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where AS 


Id 


Is the Ideal entropy of mixing of the 
X , and AS„u Is the Ideal entropy of mixing of the 

a |3 n 

solution at composition X^. 
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solid solution at 

B solid 


LIF-HgFg . - The L1F-MgF2 phase diagram (fig. 2) demonstrates extensive 
solid solubility for both terminal solid solutions; however the exact limits 
are not clear at the eutectic temperature. For our calculations the following 
values were used: the L1F rich solid solution has 77 percent L1F,1 and the 

MgF 2 rich solution has 61 percent L1F. The heat of mixing for the L1F-MgF2 
melt has been determined by Hong and Kleppa (ref. 8) at 1354 K, and It only 
showed a slight variation with composition: 
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AH M 
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MgF, 


-24.09 - 2.36 X, 


MgF, 


1.025 X, 


MgF, 


kJ/mol 


(26) 

No data exists as to the heat of mixing In the solid solutions; however they 
can be estimated by combination of phase diagram Information with the heat of 
mixing data for the liquid and the Gibb's free energy of fusion of the compo- 
nents at the eutectic temperature. Assuming Ideal entropies of mixing for 
both the solid and liquid solutions, equation (9) becomes 


o(T e ) 


RT an X* +■ RT an y* - RT an a^ = -AG f 1 


(27) 


where the activity coefficients of L1F and MgF 5 In the melt can be derived from 
equations (11) to (13) and (26) as c 


RT an = -24091 X 2 


MgF, 


2363 


(' - 2 X L1f) X MgF, 


1025 


(l - 3 X UF )X 


RT ln yLc = -24091 X 2 - 4727 X 2 , • X M „ 

MgF 2 MgF 2 LI F MgF^ 


3075 X 


L1F 


LI F/ MgF 2 
( 28 ) 

< 2 

MgF 0 


(29) 
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All compositions are given In mole percent. 
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Based on the equivalent Ionic fraction approach the activities of LI F and MgFp 
In the solid solution are 



(30) 

(31) 


with X 
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' ' being the equivalent Ionic fractions of L1F and MgF 0 , 
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respectively In the solid solution, l.e. 
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The activity coefficient of Li F In the solid solution In the L1F side of the 
phase diagram can be computed from equations (27), (28), and (30), and the 
activity coefficient of MgF 2 In the solid solution In the MgF 2 side of the 
phase diagram can be calculated from equations (27), (29), and (31). Knowl- 
edge of the activity coefficient of one component then can be used to deter- 
mine the heat of mixing In the solid solutions with the assumption of a regular 
solution model (ref. 11). In this case the higher order terms In equation (5) 
are neglected and the heat of mixing per equivalent of the solution Is 


AH M 


(S) 



(34) 


For the regular solution approximation the activity coefficient Is 


RT 8,n y 


(S) 

1 



and the heat of mixing/mole of solid solution Is 


- (’ ♦ W 2 ) • - M ‘ s > 


(35) 


(36) 


From a knowledge of the activity coefficient of a component In a solution the 
constant "A 1 " can be calculated; then the heat of mixing In the solid solution 
can be obtained from equations (34) and (36). With values for the heat of 
mixing In the liquid and the solid solutions, equation (24) can be utilized, 
and for L1F-MgF 2 eutectic at 1013 K the heat of fusion via this approach Is 
20.878 kJ/mole. 
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With the assumption of random mixing In both the liquid and solid phases 
and using the equivalent Ionic fractions for the solid solutions, the entropy 
of fusion of the L1F-MgF2 eutectic was calculated to be 26.5 J mol -1 which 
leads to a heat of fusion of 26.8 kJ/mole. Comparison of this estimate to 
that calculated from the heat of mixing considerations shows a large discrep- 
ancy (~25 percent) between the two methods which arises from the many assump- 
tions necessary to compute the terms In equations (1) and (2). 


KHgF 3 ~KCaF 3 . - This system exhibits an eutectic equilibrium among liquid 
containing ~60 percent KCaF 3 and solid phases with 20 and 90 percent KCaF 3 at 
1258 K (fig. 4). Since the heat capacity of the liquid and the constants In 
equation (4) for heat capacities of the solids are not known, the heat of 
fusion for KMgG 3 and KCaF 3 will be assumed to be Independent of temperature. 
Similarly there exists no data for the heats of mixing for either the liquid 
or solid solutions; however they can be estimated from the phase diagram 
according to equations (9) to (12) and the assumptions that: (a) random mix- 

ing exists In all solutions, (b) the regular solution approximation holds for 
the melt, and (c) In the solid solutions Raoult's law Is followed by the sol- 
vent and Henry's law Is applicable to the solute. For the equilibrium between 
the eutectic melt and the KCaF 3 rich solid solution (8) equation (9) becomes 


RT an X 


e(D 

KCAFg 



A' - RT an X 
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= -AG 


o(T e ) 
f ,KCaF, 


(37) 


The solution of equation (37) yields A' to be -556 J/mol . Thus the heat of 
mixing for the melt becomes 



-556 X 


a 

KCaF 3 



J/mol 


(38) 


From Raoult's law the activity of KMgF 3 In the KMgF 3 rich solid solution (a) 

Is ; hence the thermodynamic equilibrium between the solid solutions 

3 

a and 8 requires 


*KMgFj * X KH 9 F 3 * 4gF 3 < 39 > 

Solution of equation (39) gives the activity coefficient of KMgF 3 In the KCaF 3 
rich solid solution (8). Likewise the activity coefficient of KCaF 3 In the 
a solid solution can be determined. With a random distribution of cations, 
the heat of mixing/mole of each solid solution can be expressed as 


4H S • X KCaF 3 
4H S ' X KHgF 3 


RT an 
RT an 


a 


Y KCaF 3 

8 

Y KMgF 3 


(40) 

(41) 


which when evaluated give 3.14 and 2.17 kO/mol for a and 8 solid solutions, 
respectively. Substitution of the appropriate terms Into equation (24) yields 
a heat of fusion for the KCaF 3 -KMgF 3 eutectic of 58.6 kJ/mol. This can be 
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contrasted with a value of 65.7 kJ/mol which was determined from equation (25) 
assuming random mixing In all solutions. 

NaHgF^ - KMgF^ . - While NaMgF 3 and KMgF 3 are completely soluble In one 
another (fig. 4) for both liquid and solid phases, the two phase equilibrium 
between liquid and solid possesses a minimum at 30 percent KMgF 3 and 1288 K 
where congruent melting takes place. Although no data exist for the heats 
of mixing In either the liquid or solid solutions, It Is expected that 

AH 8 > AH* > 0 (ref. 17). An estimate of the heat of mixing for the solid 

solution can be obtained from the phase diagram and the assumption of Ideal 

0 

mixing for the melt: l.e., AH h = 0. With the further assumption that (a) the 

heat of fusion Is Independent of temperature, due to the lack of Information on 
heat capacities, and (b) random mixing exists In the solid solution, 
equation (9) becomes 


<Sf, K) - 1288 R • T| 


NaMgF, 


. r o(l 288 K) , OQQ D . S 
AG f,KMgF, = 1288 R * *" y KI 


KMgF, 


Thus the heat of mixing of the solid solution at X 


KMgF, 


0.3 Is 


S(0.3 KMgF ) - 

ah m = 901 R an T NaHg p 3 + 386 R an T K ^ gF ^ 


(42) 

(43) 


(44) 


Combinations of equations (42) to (44) gives the heat of mixing for the solid 
solution at 30 mol percent KMgF 3 to be 1.4 kJ/mol, and this value can be sub- 
stituted Into the expression for the enthalpy of fusion 


AH 


0.3 KMgF 3 
f 


0.7 AH 


1288 K 
f, NaMgF 
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+ 0.3 AH 


1288 K 
f.KMgF 


3 


S( 0 . 3 KMgFJ 

4H « 


(45) 


which yields a heat of fusion at the congruent melting point to be 70 kJ/mol . 

Because the solid and liquid phases are of the same composition at the 
congruent melting point, the entropy of mixing would be the same for both 
liquid and solid solutions If random mixing Is assumed. Therefore the entropy 
of fusion at the congruent melting point Is given by 

0.3 KMgF 

AS f = 0,7 AS f,NaMgF 3 + 0,3 AS f,KMgF 3 (46) 

Evaluation of this equation yields 54.3 J K _1 mol -1 which converts to a heat 
of fusion of 70 kJ/mol from equation (2a). 
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Binary Eutectic Systems With Congruent Intermediate Compounds 


Theory . - Figure 6 Illustrates a binary system with one congruent melting 
compound at composition X c and two eutectics at X e -| and X e 2 - Blrchenall 
and Rlechman (ref. 13) have outlined a method to estimate the heat of fusion 
for this type of system where there Is negligible solid solubility of one com- 
ponent In the other. For the A-rlch eutectic at X e -| , 1 mol of liquid upon 
freezing Is converted to X e -j mol of compound c and (1 - X e -|/X e ) moles of 
pure A. Similarly for the B-rlch eutectic at X e 2 , 1 mol of liquid upon 
freezing Is converted to (1 - X e 2 ) moles of the Intermediate compound e and 
1 - {(1 - X e 2 )/X c } moles of pure B. The heat of fuslon/mole of the eutectic 
Is given by 

■ d - x ei> 4H M * x el sh M * sh H 1(1> - x el 4H ! <« 7 > 

T T 

ah' 2 - (1 - X e2 > AH f ' 2 ♦ K e2 4H f ' 2 . AH' 2 < 1 > - (1 - X e2 > AH* (48) 

where AH® 1 ^ and AH® 2 ^ are the heat of mixing of the eutectic liquid 
at compositions el and e2, respectively; AH^ Is the heat of formation for 
1 mole of the Intermediate compound e from the components A and B. 

Disregarding the entropy of mixing for the Intermediate compound c, the 
entropy of fuslon/mole of the eutectic can be calculated In the same manner as 
binary eutectics without terminal solid solutions (eq. (17)). In general sig- 
nificant negative deviation would be expected In the liquid state for systems 
which form Intermediate compound solid phases; thus the actual entropy of mix- 
ing should be lower than that calculated for an Ideal solution. Because of 
this, the excess entropy of mixing must be taken Into account; fortunately for 
systems where the heat of mixing Is known, the excess entropy of mixing can be 
estimated from phase diagrams via equations (9) to (14). 

N aF-MqF 2 . - The compound NaMgF 3 Is formed In the NaF-MgF 2 system (ref. 18) 
with eutectics forming on the NaF side of the phase diagram at 23 percent MgF 2 
and 1103 K and on the MgF 2 side at 60 percent MgF 2 and 1273 K. Hong and Kleppa 
(ref. 8) have measured the heat of mixing of the NaF-MgF melt at 1345 K for 

0 < v 2 < °- 7 - 


NaF-MgF 

AH M 


= X NaF 



^-48.58 - 20.309 



+■ 8.732 X 


2 

MgF 



kJ/mol 


(49) 


They also estimated the enthalpy of formation of NaMgF 3 to be -12.13 kJ/mol 
at 1300 K. Assuming that this value Is unchanged by temperature, the heats of 
fusion of the two eutectics calculated from equations (47) and (48) are 
29.8 kJ/mol of NaF-23MgF 2 and 38 kJ/mol of NaF-60MgF 2 . 
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The experimental llquldus on the NaF side of the phase diagram and the 
eutectic point (60 percent MgF 2 , 1273 K) on the MgF 2 side were used with 
the heat of mixing data (eq. (49)) to calculate the excess entropy of mixing 
via equations (9) to (14). from regression analysis 


: E(i) 


= X 


NaF 


MgF, 


(' 


18.99 - 47.23 X 


MgF, 


+ 43.59 X 


MgF. 


JK 1 mol 


-1 


(50) 

Combination of equation (50) with equation (17) to determine the entropies of 
fusion for the eutectic yields 28.4 JK -1 mol~l and 31.15 JK-1 mol-1 for the 
NaF- and MgF2-r1ch eutectics, respectively. Substitution of these values Into 
equation (2a) gives the heats of fusion of 31.3 arid 39.7 kJ/mol for the NaF- 
and MgF2~r1ch eutectics, respectively. These values are In close agreement 
(~5 percent) with those calculated directly from the heat of mixing data. 


KF-MgF 2 . - The Intermediate phase KMgF 3 forms an eutectic with KF at 
1051 K at 14 percent MgF 2 and an eutectic with MgF 2 at 1281 K at 68.7 percent 
MgF 2 (ref. 19). The heat of mixing for KF-MgF 2 melt was found to be (ref. 8) 


KF-MgF ? ( 8.) 

AH M = X KF 


-63.082 - 64.751 X„ c + 37.75 X‘ 

l*Qr ^ nQr ^ MQr ^ 


kJ/mol 

(51) 


and the enthalpy of formation of KMgF 3 from KF and MgF 2 was estimated to. 
be -22.6 kJ/mol at 1130 K. Assuming that the heat of formation for this 
compound Is Independent of temperature, as no data for the heat capacity of 
KMgF 3 exists; calculations based on equations (47) and (48) give 25.9 and 
35.1 kJ/mol for the heats of fusion for the KF- and MgF2-r1ch eutectics, 
respectively. 


Estimation of the excess entropy of mixing along the llquldus curves for 
the KF-rlch side of the phase diagram gave positive values which seem unrea- 
sonable In light of the strong Interaction between KF and MgF 2 , as evidenced 
by the large negative values for the heat of mixing. Because Hongand Kleppa's 

measurements ceased at X u c = 0.7, the required extrapolations to higher 

MgFp 

MgF 2 concentrations could lead to error In the entropy of mixing values; 
thus the excess entropy of mixing was calculated using only two points (50 and 
68.7 mol percent MgF 2 ). From the literature (ref. 8) the entropy of fusion 
for KMgF 3 Is 48.53 JK -1 mol -1 which Is related by 

1d( 0 . 5 MgF ) E(0.5 MgF ) 

AS f,KMgF 3 = AS f,KF * AS f,MgF 2 + 2 AS M * 2 S M (52) 


to the entropies for KF and MgF^ (A ^ and 

0.5 MgF (Id) 

Ideal entropies of mixing, AS^ and 

composition. After solving equation (52) for 


AS 


f ,MgF 


) plus the excess and 


E(0. 5 MgF ) 

S M , at the equimolar 

E ( 0 . 5 MgF ) 

S M and combining this 
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expression with the excess entropy of mixing for 68.7 mol percent MgF 2 from 
equations (9) to (14), the excess entropy of mixing for the melt Is 


cEU) 




^-25.85 


5.52 X 


MgF, 


JK -1 mol" 


( 53 ) 


From this equation and the entropies of fusion for KF and MgF 2 , equation (17) 
yields 23.6 and 27.5 JK _1 mol -1 for the entropies of fusion of the KF- and 
HgF 2 ~r 1 ch eutectics, respectively. In turn these values give (eq. (2a)) 
enthalpies of fusion of 24.8 and 35.2 kJ/mol for the respective eutectics 
which are In very good agreement with the values calculated directly from the 
heats of mixing. 


KF-CaFg . - The KF-CaF 2 system (ref. 20) has an Intermediate compound at 
the equimolar composition which melts at 1343 K and possesses two eutectics: 
at 1055 K with 15 percent CaF 2 and at 1333 K with 62 percent CaF 2 . Kleppa and 
Hong (ref. 6 ) have measured the heat of mixing for this system where 


KF-CaFp( 9.) 

AH m = X kf • X CaF (-33.848 + 1.799 X CaF \ kJ/mol (54) 


Neither the heats of formation nor fusion are known for KCaF 3 *, however these 
values can be estimated from the excess entropy of mixing. From the llquldus 
curves on both sides of the phase diagram, equations (9) to (14), and the heat 
of mixing expression (eq. (54)), the excess entropy of mixing In the melt Is 


cE(D 



X CaF 2 (-13.8*1.51 




mol 


(55) 


The entropy of fusion for KCaF 3 can be obtained from an expression modeled 
after equation (52) and Is computed to be 50.85 JK _1 mol -1 ; thus the enthalpy 
of fusion for this Intermediate compound Is 68.3 kJ/mol . The heat of forma- 

f 

tlon of KCaF 3 from KF and CaF 2> AH KCa(r , at Its melting point can be calculated 
from 


AH 


0.5 CaF (H) 

KCaF 3 = AH f,KF * &H f,CaF 2 + 2AH M " AH f,KCaF 3 


(56) 


and gives -19.53 kJ/mol . 

The estimated heats of fusion for the two KF-CaF 2 eutectics were Inde- 
pendent of the method of calculation and were obtained to be 27.6 kJ/mol for 
KF-l 5Caf 2 and 33.2 kJ/mol for KF-62CaF 2 . 

MF-AIF^ . - Three metal fluoride (M = Na, LI, or K) -AIF 3 phase diagrams 
have an Intermediate compound M 3 A 1 F 5 at 25 percent AIF 3 (refs. 21 to 23) which 
forms an eutectic with pure MF. As these systems have not been extensively 
studied beyond the f^AlFg composition, only the MF-rlch eutectics can be 
considered for energy storage applications. Because of Its volatility, the 
melting point and heat of fusion for AIF 3 are not known; however Hong and 
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Kleppa (ref. 7) have estimated the heat of fusion to be 112 kJ/mol at 1298 K. 
Since the heat capacity data for liquid MF 3 are not available, the enthalpy 
of fusion for this compound will be assumed to be Independent of temperature. 
Heats of mixing for all three MF-A 1 F 3 systems have been measured by Hong and 
Kleppa (ref. 7) and their results are presented In figure 7. This data com- 
bined with the enthalpies of formation of M 3 AIF 5 , as computed from the thermo- 
dynamic data In Barln and Knack's text (ref. 4), leads to the following heats 
of fusion for the three eutectics: 25.8 kJ/mol for NaF-1 4A1 F 3 at 1161 K, 

28 kJ/mol for KF-6.8A1F 3 at 1113 K, and 24.7 kJ/mol for L1F-14.5A1F 3 at 983 K. 
Since the entropy of fusion of A1 F 3 Is not known, no attempt was made to cal- 
culate the heat of fusion via equations (17) and (2a). 


TERNARY EUTECTIC SYSTEMS 


The ternary system NaF-CaF 2 MgF 2 exhibits two ternary eutectics (ref. 24) 
at NaF-27.2CaF2-36.5MgF 2 and 1178 K and at NaF- 22 . 7 CaF 2 -l 2 . 8 MgF 2 and 1018 K 
which can be analyzed via equation (15). The heat of mixing for the NaF-CaF 2 
and NaF-MgF 2 can be calculated from equations (23) and (49), respectively. No 
known data exists for the heat of mixing for CaF 2 -MgF 2 ; however with the 
assumption of the regular solution model, the heat of mixing can be calculated 
from the binary phase diagrams where 


CaF_-MgF ( 1 ) 
AH M 


-11.26 




kJ/mol 


(57) 


Substitution of the binary heat of mixing values along constant composition 
paths calculated from equations (23), (49), and (57) In equation (15) gives 

NaF-CaF -MgF (l) 

the heat of mixing AH eM for the ternary eutectic melt. Use of 

NaF-CaF -MgF ( 8 . ) 

AH gM along with the heat of fusion of the Individual components 

at the eutectic temperature In 


AH ef - X NaF aH flNaF + X CaF, 


fiHr p r * X M r 
f,CaF 2 MgF 2 


AH 


f.MgF 


2 


+ AH 


NaF-CaF 2 -MgF 2 (!l) 

eM 


(58) 


yields the heat of fusion for the ternary eutectics to be 30.75 and 

26.94 kJ/mol for NaF-27 . 2CaF 2 -36 . 5MgF 2 and NaF-22.7CaF 2 -12.8MgF 2 , respectively. 


CONCLUDING REMARKS 

In space power applications the heat of fusion per unit mass Is probably 
more Important than heat/mol since the energy storage media must be put Into 
earth orbit. On the basis of literature data and thermodynamic calculations 
there appear to be a large number of fluoride systems (table III), Including 
pure compounds and eutectic mixtures, which have relatively high enthalpies of 
fusion and melting points between 980 and 1378 K. Unfortunately no suitable 
systems have been Identified between 1378 and the 1536 K melting point of MgF 2 - 
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TABLE I. - THERMODYNAMIC DATA FOR PURE SALTS 


Compound 


Enthalpy of 
fusion, 
kJ/mol 

Entropy of fusion, 
JK" 1 mol" 1 

The constants for 
capacity of solid ( 
JK -1 mol" 

the heat 

jq. (4)), 

Heat capacity 
of the melt, 
OK"-*- mol"-*- 

Reference 

A 

B 

C 

a CaF 2 

1691 

29.706 

17.568 

107.98 

10.46 







(1424 to 1691 K) 

(1424 to 1691 K) 











99.997 

4 




20.92 

59.83 

30.459 

1.966 






(<1424 K) 

(298 to 1424 K) 



KF 

1130 

28.242 

24.99 

45.98 

14.42 


66.944 

4 

LiF 

1121 

27.196 

24.254 

42.689 

17.417 

-5.3 

64.18 

4 

LaF 3 

1766 

50.245 

28.45 

77.73 

19.953 


152.716 

4 

MgF 2 

1536 

58.157 

37.86 

70.83 

10.543 

-9.2 

94.43 

4 

NaF 

1269 

33.137 

26.112 

46.588 

9.899 

-2.129 

70.567 

4 

b YF 3 

1428 

27.97 

19.585 

-319.44 

212.798 

2817.92 






(1350 to 1428 K) 

(1350 to 1428 K) 



1 








133.662 

4 




43.60 

99.411 

7.426 

-5.69 






(<1350 K) 

(<1350 K) 




KMgF 3 

1343 

65.178 

48.53 





8 

NaMgF 3 

1303 

74.14 

56.9 




— 

8 

c KCaF 3 

1343 

68.299 

50.85 





- 


a Phase change in solid at 1424 K, AHt rans = 4.769 kJ/mol. 

b Phase change in solid at 1350 K, AH^ rans = 32.426 kJ/mol. 

c Enthalpy and entropy of fusion estimated from the present calculations. 














TABLE II. - ESTIMATED HEATS OF FUSION FOR THE SIMPLE BINARY EUTECTIC SYSTEMS 


[Compositions are In mole percent, e.g., A - mol percent B. NA - data 
on heat of mixing not available.] 


Composition, 
mole percent 

Eutectic 

temperature, 

K 

Heat of fusion 
estimated from 
equation (16), 
kJ/mol 

Heat of fusion 
estimated from 
equation (17), 
kJ/mol 

Heat of fusion 
estimated from 
equation (20), 
kJ/mol 

KF-40NaF 

994 

28.233 

29.2 

29.15 

Lif-22CaF 2 

1039 

27.86 

28.17 

28.47 

NaF-32CaF 2 

1083 

28.04 

30.95 

30.35 

MgF 2 -70YF 3 

1238 

NA 

58.5 

54.8 

CaF 2 -50MgF 2 

1253 

NA 

44 

45.9 

MgF 2 -37LaF 3 

1378 

NA 

54.9 

48.5 









TABLE III. - HEAT OF FUSION PER UNIT HASS FOR SEVERAL 

FLUORIDE SALTS 


Composition, 
mole percent 

Temperature, 

K 

Heat of fusion, 
kJ/gm 

LI F-l 4 . 5A1 F 3 

983 

0.721 

NaF-60KF 

994 

.546 to 0.564 

L1F-33MgF 2 

1013 

.55 to 0.7 

NaF-22CaF 2 -13 MgF 2 

1018 

.51 

L1F-22CaF 2 

1039 

.745 to 0.761 

KF-14HgF 2 

1043 

.42 to 0.44 

KF-1 5CaF 2 

1055 

.452 

L1 3 A1F 6 

1059 

.54 

NaF-32CaF 2 

1083 

.52 to 0.56 

NaF-23MgF 2 

1103 

.64 to 0.67 

KF-6.8A1F 3 

1113 

.47 

LI F 

1118 

a l .087 

KF 

1129 

a .486 

NaF-1 4A1 F 3 

1161 

.539 

NaF-27CaF 2 -36 MgF 2 

1178 

.52 

MgF 2 -70YF 3 

1238 

.45 to 0.48 

CaF 2 -50MgF 2 

1253 

.61 to 0.65 

KCaF 3 -40KHgF 3 

1258 

.45 to 0.5 

NaMgF 3 -30KMgF 3 

1259 

.641 

K 3 A1F 6 

1263 

a .47 

NaF 

1268 

a .789 

NaF-60MgF 2 

1273 

.7 to 0.732 

KF-68.7MgF 2 

1281 

.57 to 0.58 

Na 3 AlF 6 

1283 

a .53 

NaMgF 3 

1303 

b .71 1 

KF-62CaF 2 

1333 

.47 

KMgF 3 

1343 

b . 541 

KCaF 3 

1343 

.501 

FeF 2 

1373 

a . 553 

MgF 2 -37LaF 3 

1378 

.43 to 0.49 

MgF 2 

1536 

a .933 


a Obta1ned from reference 4. 
^Obtained from reference 8. 
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Figure 1. - Definition of terms used to calculate the enthalpy of mixing for 
ternary liquid solutions from binary heat of mixing data. 



Figure 2. - Schematic of a phase diagram having a single eutectic 
with terminal solid solutions. 




NaMgFj 


Content, mol % 

Figure 5. - NaMgFj - KMgFj phase diagram (ref. 16). 









Figure 6. - Schematic of a phase diagram with a congruently 
melting intermediate compound and two eutectics. 
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